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In this paper we propose and analyze a fully discrete method for a direct bound- 
■ ary integral formulation of the scattering of a transient acoustic wave by a sound-soft 

obstable. The method uses Galerkin-BEM in the space variables and three different 
choices of time-stepping strategies based on Convolution Quadrature. The numeri- 
cal analysis of the method is carried out directly in the time domain, not reverting 
to Laplace transform techniques. 
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§ ■ 1 Introduction 

In this paper we propose and analyze a fully discrete method for the direct boundary in- 
tegral formulation of the Dirichlet problem for the causal acoustic wave equation, exterior 
to a domain with Lipschitz boundary in M. d (d = 2 or 3). The method arises from using 
a general Galerkin semidiscretization-in-space and multistep-based Convolution Quadra- 
ture (CQ) in time. From the point of view of the numerical method, this paper extends 
work in [12] and [3]. A survey of recent results for CQ-BEM discretization of a wide 
variety transient problems can be found in [TJ. 

Analytical literature on time-domain integral equations has typically been focused 
on integral equations of the first kind arising from indirect formulations. The origin of 
these techniques was based on Galerkin- in-time methods [1], while CQ techniques were 
developed only about one decade later. The present paper uses a direct formulation, 
leading to an integral equation of the first kind similar to those treated in [TJ [121 E] • The 
main differences lie in the fact that data appear under the action of a retarded integral 
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operator (which will have to be discretized as well) and that the unknown on the boundary 
is a quantity of physical interest. 

In this paper we propose the development of a systematic analysis of CQ-BEM taking 
care of all aspects of discretization: (a) data interpolation on the computational grid, (b) 
Galerkin semidiscretization-in-space of the associated retarded integral equation, (c) dis- 
cretization in time (using CQ) of the integral operators in both sides of the equation, (d) 
discretization in time of the postprocessed potentials (acting on the data and unknown 
of the integral equation) to obtain the scattered wave field in exterior points. The main 
difference with the traditional black-box analysis proposed by Lubich [T2] is in the fact 
that we propose to do most of the analysis directly in the time domain. Original work in 
the analysis of CQ-BEM dealt only with the simplest retarded boundary integral equa- 
tions, that are coercive in the Laplace domain. Coercivity is inherited by the Galerkin 
semidiscrete-in-space problem, but some properties of the fully discrete problem (includ- 
ing postprocessing of the solution to obtain the associated potentials and treatment of 
data that appear under the action of retarded integral operators) have to be investigated 
in a more direct fashion [10] . More recently, some estimates in the time-domain [8j [16] 
have expanded the analytical toolbox that can be used to prove error estimates for full 
discretization of retarded boundary integral equations. It has to be noted that most of 
the literature that is relevant for this analysis had been carried out using Laplace trans- 
forms -the paper [15] seems to be a lone exception-. The passage through the Laplace 
domain makes for a relatively streamlined analysis that can be expanded to a wide va- 
riety of problems [TOj [6] but is likely to yield less sharp results than a direct analysis in 
the time domain. As announced, in this paper we will develop the more recent technol- 
ogy of time-domain estimates to show properties of the Galerkin semidiscrete-in-space 
problem (these are pertinent for other kind of time- discretization methods, using Galerkin 
schemes) and of the full discretization of the problem. In particular, we will obtain a 
proof of convergence of the trapezoidal rule CQ method that is not directly reachable 
in the Laplace domain. The tools for this analysis are varied but not complicated: (a) 
identification of the weak convolutional retarded integral equations and layer potentials 
with strong solutions of problems in finite domains for finite time intervals, (b) interpreta- 
tion of Galerkin semidiscretization-in-space with exotic transmission problems following 
[T0| IT6] , (c) use of the well understood theory of Co-groups of isometries [H] to obtain 
estimates for the resulting dynamical systems [8], (d) understanding of the process of CQ 
time-discretization as a direct discretization of the exotic transmission problems in the 
time-domain (the essence of this idea is already present but not exploted in [TJ]) and (e) 
application of standard techniques for numerical analysis of the wave equation in bounded 
domains to work out the analysis of the fully discrete method . 

Foreword. Elementary properties of basic Sobolev spaces if x (f2) and i? ±1,/2 (r), the 
trace operator and the weak normal derivative, will be used without further reference. The 
pertinent results can be found in any advanced textbook of elliptic PDE. The monograph 
[T3] contains all of them, as well as some results about steady-state layer potentials and 
integral operators that will be similar to the ones we will be developing in this paper 
(and that are used to prove background results that will be explicitly mentioned as they 
are used). While possible, we will make an effort in clarifying the source of constants in 
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estimates. Once this is not practical any more, we will use the convention of admitting 
C > to be a constant indendendent of the associated discretization parameters (h and 
k in this paper). Vector- valued distributions appear in the background of the theory of 
retarded layer potentials and integral operators. Their use has been outsourced to some 
preliminary papers PUJ El [TB] that relate strong and weak solutions of the wave equation. 
Here we will only employ the basic idea of a causal distribution with values on a space 
X as a sequentially bounded map X?(R) — > X that vanishes when applied to elements of 
T>((— oo,0). The concepts of differentiation and Laplace transform are then identical to 
those of scalar distributions. 

2 An integral formulation of the scattering problem 

Let Q~ be a bounded open set in M. d (with d = 2 or 3) with Lipschitz boundary T. 
We admit the possibility that f2~ is not connected, but we demand the complementary 
domain Q + := M. d \ Q to be connected. The problem of scattering of an incident wave 
by a sound-soft obstacle can be written by means of the Initial Boundary Value Problem 

u tt = Am in Vt > 0, (2.1a) 
u + u inc = onT, Vt>0, (2.1b) 
u{ ■ ,0) = u t ( ■ ,0) = on T. (2.1c) 

The incident wave u mc is a known function. For the model equation to be meaningful we 
have to assume that u mc ( ■ , 0) = uf LC { ■ , 0) = in a neighborhood of V . The unknown 
in (12.1 p is the scattered wave field, while the total wave is u + u mc . There is no need to 
impose a radiation condition at infinity since causality of the wave equation takes care of 
the fact that the support of the solution of (I2.ip . for any given t, is compact. 

For all purposes (expository and analytic), it is convenient to understand functions of 
the space and time variables as functions of t with values on a space of function. Therefore, 
instead of considering u = w(x, t) as a function in Q + x [0, 00), we will consider u = u(t), 
where u(t) G H 1 (Q + ) for all t. It will also be convenient to refer to causal functions as 
functions £ : R — >■ X (here X is any Hilbert space) such that = for all t < 0. The 
concept of causality can be easily extended to distributions with values in the space X. 

An integral representation of the solution of (12. ip starts by taking the value of the 
incident wave on T for positive values of the time variable. If ^ : if 1 (fi ± ) — > iJ 1 / 2 (r) are 
the trace operators on T, we consider the causal function 

ip : R -> H 1/2 (T) such that (p(t) = 7 + w inc ( • , t) Vt > 0. (2.2) 

Note that the required regularity of the incident wave for this process to be meaningful is 
local H 1 behavior in a neighborhood of T and that u mc can have singularities away from 
the scattering boundary (this is the case for waves originated by acoustic sources). 

Consider now the single and double layer retarded acoustic potentials. Their strong 
expressions in the three dimensional case (valid for smooth-in-space densities written as 
functions of the space and time variables) are 

(.S*A)(x,t) := / ^V^^^ ^y) xGR 3 \T, (2.3) 
J r 47r|x-y| 
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and 



■Ky)dr(y), xel d \r, (2.4) 

z=y 

respectively. In (12. 4 ft the vector u(y) denotes the unit outward pointing normal vector 
at the point y G T. The notation for the layer potentials in H 2 . 3 f) and (12 Ah uses the 
convolutional symbol to emphasize the fact that these are time-convolution operators 
(see [9] for a rigorous introduction of these operators in the sense of distributions), since 
we will take advantage of this convolutional structure for the discretization in the time 
variable. 

By using direct arguments in the time domain [5] or employing Laplace transforms 
[TJ |2], it is possible to prove that if A is a causal distribution with values in the space 
iJ -1 / 2 (r), then S * A is a causal distribution with values in the space 

#i(R d \ T) := {v G H\R d \ T) : Am G L 2 (M. d \ T)}, (2.5) 

satisfying 

A(S * A) = ^(<S * A) inR d \r (2.6) 

and 

[ 7 (<S*A)] := 7 + (S*A)-7 + 0P*A) = 0, [d v (S*X)] := <9+(<S*A) -(9~(<S*A) = A. (2.7) 

Note that the Laplace operator (in the sense of distributions in R^ \ T) and the exterior 
and interior normal derivatives are well defined in the space H^(M. d \T). Expressions (12.61) 
and (12. 7p can be understood as equalities of causal distributions with values in L 2 (R d \r), 
H l l 2 {T) and H~ l l 2 (T) respectively. The second derivative in (I2.6P is defined in the sense 
of vector valued distributions. The first of the jump relations (12. 7p allows us to define the 
retarded integral operator 

V* A :=7 + (S* A) =7"(<S* A), (2.8) 

whose integral expression in the three dimensional case (for smooth enough densities) 
coincides with that of S * A, with x G V now. 

If (p is a causal distribution with values in if 1/,2 (r), then V * is a causal distribution 
with values in i?^(R d \ T) satisfying 

A(V*cp) = £z(V*ip) inR d \r (2.9) 

and 

[ 7 (Z> * ?)] = [a,(2>*y01=O. (2.10) 
We then define the retarded boundary integral operator 

£*£ = I 7 +(P*£) + | 7 -(D*£). (2.11) 

An integral expression for this operator in the three dimensional case coincides with that 
of the layer operator (see (I2.4p ). Any causal if^(R d \r)-valued tempered distribution 
u such that 

u = Au in R d \ T 



(P*^)(x,t) 



(f{z,t- \x-y\) 
47rlx — vl 
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(with equality as L 2 (IR d \ r)-valued distributions and with the usual notation u = 4^r) 
can be represented with Kirchhoff 's formula 

u = S * \d v u\ — V * [77/]. 

Therefore, if we consider a solution of (12. ip as a causal tempered if 1 (f2 + )-valued distri- 
bution that is extended by zero to Q~ and denote <p as in (12. 2p . we can write 

u = -S*\-V*ip (2.12) 

where A := d^u. Using the definitions of the boundary operators in (12.81) and (12.111) 
as well as the first of the jump relations (12.101) . it follows that A satisfies the following 
equation 

V*\ = l(p-lC*(p. (2.13) 

The analysis of [1] includes a proof of the unique solvability of the operator equation in 
(I2.13P and a Sobolev estimate for the solution of V*£ = ip. Also [SJ Theorem 6.2] contains 
an estimate of the solution operator for equation (I2.13P and its postprocessing (I2.12p . 



3 Discretization 

We start by assuming that the data function ip has been approximated. This is the 
usual approach of the engineering literature (see the exposition of a very similar family of 
methods for elastic waves in [17] for instance) and will be for us a motive to studying the 
propagation of errors in data, a study that will be needed for analysis of the fully discrete 
schemes. We therefore assume that a causal function cp^ : R — > H l l 2 (T) is given as an 
approximation to (p. 

3.1 Semidiscretization in space 

We consider a discrete space X h C L°°(T) and substitute (12.131) by the search of a causal 
function A^ : K — > such that 

(fjL h , V* X h )r = (t-ih, \<Ph - K, * ip h ) r Wfi h eX h . (3.1) 

Here and in the sequel the angled brackets denote the if~ 1/,2 (r) x H lj/2 (T) duality product. 
The solution of (13.11) is then used for the discrete representation formula 

u h := -S * X h - V * (p h . (3.2) 

For the sake of clarity, let us write down the system (13.11) in the three dimensional 
case, when data have been approximated by a function ip^ : M — > C W 1,0O (T), where 
Yh is finite dimensional. Let {N 1} . . . , Nj} and {Mi, . . . , M^} be respective bases of X^ 
and Yh. Data and unknown can then be represented by their coefficients: 

K J 

<p h (t) = <Pk(t)M k , A h (t) = X At) N r ( 3 - 3 ) 
fc=i j=i 
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This is equivalent to substituting the Yh- and Xh- valued functions by a finite set of casual 
scalar functions. Problem (13 .ip is equivalent to the system 

tf /f^A,( t -|x-y|)dr(x)dr ( y) 
j^JrJr 47r|x-y| 

+ W f M i , fc (x,y)( yfc(t ~j X ~ y|) + <p k {t- [x-y[))dr(x)dr(y) (3.4) 
k=1 Jr Jr l x y| 

(for i — 1, . . . , J), where 

M i>fe (x,y) := (X ~ y) Jy ( 2 y) iV,(x)M fc (y). 

The following spaces will be relevant in the sequel: 

W fe (M;X) := {peC^iR-X) : supppC [0,oo), g4(M;I)}, 
C fc (M;X) := {p £ C k (R; X) : supp p C [0, oo)}. 

Theorem 3.1. Let y? £ W 4 (M; if 1/2 (r)) ; assume that the solution of (I2.13P satisfies 
A £ W 2 (M;#- 1/2 (r)) and Zet u be given by fl2J2|) . ^Issnme a/so tfiat ^ £ W^M;!^). 
Then the semidiscrete equation (13. ip has a unique solution Xh- Let finally Uh be given by 
(I3~2l) . Then, for all t > 0, 

IKA-AfcXOU-w < C(l + t)(^/ t ||(AW-n fc AW)(r)||_ 1/2>r dr 

V £=0 ^ 
4 

+ Ell^ ) -^ ) )( r )llv 2 ,rdr) (3.5) 



£=0 

(«-«a)(*)IIi*< < C(l + t)(f\ /'||(AW-n A AW)(r)||_ 1/2 ^T 

V £=0 ^ 



+ Ell(^ ) -^ ) )(r)||i/ 2 ,rdr), (3.6) 

£=0 



where Uh '■ H 1//2 (r) — >• zs t/ie orthogonal projection onto Xh- 

Note that by [8, Theorem 6.2], if p £ W 4 (M; # 1/2 (r)), then A £ Cg(R; #~ 1/2 (r)). 
Time regularity of the solution is then guaranteed by the sufficient (but not necessary) 
condition <p £ W^(M; H 1 / 2 ^)). The proof of Theorem O will be given at the end of 
Section [5j If data are not discretized the second group of terms in the error estimates of 
Theorem 13.11 is not needed. 
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3.2 Full discretization 



In a final step, we substitute the four time convolutions that appear in (13.11) and (13. 2p with 
a discrete convolution based on one of the applicable Convolution Quadrature methods. 
For a fixed time-step k > 0, the CQ method applied to the discretization of ( 13. ip and 
( 13. 2 p produces (in theory) casual functions A£ : R -> X h and u£ : R -> H\(R d \T). In 
practice, these functions are evaluated in equally spaced time steps t n := nn and only 
these values of the functions are obtained. To obtain values at other times, the method 
has to be run again, starting at to := — £ t for instance. Therefore, even if the theory of 
CQ deals with functions of continuous time, in practice the solutions can be understood 
as functions of discrete time, i.e., sequences. 

Let us briefly explain what the CQ discretization of (13 ,4p consists of. First of all, 
we consider the complex matrix valued functions V^(s) G R JxJ and K/^s) G R JxX with 
elements 

Y h (s) hJ := / 1 l^yf e- |x - y| dr(x)dr(y) (3.7) 

: = / /^, fc (x,y)(l + S |x-y|)^4dr(x)dr(y) (3.8) 
Jr Jr l x ~~ y| 

and the matrix with elements 

\ k := [ iV l (x)M fc (x)dr(x). 
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Note that Yh{s) and K/^s) are the Laplace transforms of the operators that appear in 
(13 .4p . We then construct the Taylor expansions 

oo oo 

V^K-^CO) = E V ^NC m^iO) = E K £NC (3.9) 

n=0 n=0 

where 5 is one of the following functions 

{1 — Cj (backward Euler method) 

|-2C + |C 2 (BDF2) 
2 (trapezoidal rule) . 

Data discretization consists of the construction of vectors <p£[n] := (ipk(n k))^ =1 for n > 
(recall ( 13. 3ft ). The unknowns are vectors AJJ[n] G R J satisfying 

n n 

E V ^H^[n - m] = |W£[n] + E KWrtfr ~ H n > 0. (3.10) 

m=0 m=0 

We can thus associate 

j 

Afln] = (Axfn], . . . , Aj[n]) AJJ[n] := E G X, 

j'=i 
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to obtain a fully discrete approximation of \h(t n ) X(t n ). 

The postprocessing step to compute the approximated scattered field can be explained 
in a similar way. The s-domain semidiscrete single and double layer potentials correspond 
to vector valued functions with domain R 3 \ T 

B h {s) h := / ( ' ~ y) '^\ l + s\ ■ -y|)f^— ^M fc (y)dT(y). 
Jr 4vr| • -y| 2 | • -y| 

The CQ method uses the same strategy as in (13.91) to produce sequences of vector valued 
functions §£[n] and defined in R 3 \T, and uses them to construct the approximations 

n n 
m=0 m=0 

In the two-dimensional case the expressions for the fully discrete method are very similar, 
using Hankel functions instead of exponential expressions. For instance, 

Vfc(*) W := ^^^ 1) (^|x-y|)iV J (x)iV J ( y )dr(x)dr(y), 

KfcOOi,* := -7/ / ^M* - y|) (X ~ y) ' ^ (y) A^x) M fc (y) dT(x)dr(y). 
4irir l x ~y| 

Remark 3.1. When data are not approximated, the indices k have to be ignored. Instead 
of having a matrix Kft(s) we have operators i/ 1 / 2 (r) — > R J obtained by substituting the 
basis function by a general element of H l l 2 (T) . In this way, the matrix-vector products 
K£[m]<^[n — m] have to be substituted by the action of operators K K [m] : H l l 2 {T) — > R J 
on elements ip K [n — m] = ip(t n - m ) E H l l 2 {T). Similar changes have to be applied to 
the double layer potential and to the matrix \ which is now substituted by an operator 
H l l 2 (T) — > R J corresponding to testing a function with the basis functions iVj. 

The analysis of the difference between the semidiscrete and the fully discrete solution at 
the different time-steps is carried out separately for the three time-discretization methods. 
This is done in Section [61 

4 The semidiscrete Galerkin projection 

Consider a casual smooth function A : R — > H~ l l 2 (T) and a finite dimensional space 
Xh C H~ l l 2 (T). The aim of this section is the analysis of the semidiscrete discretization 
process looking for a causal function : R — >■ Xh such that 

(ji h , (V*A?)(*)) r = W,(V*A)(*)) r V^Gl,,, Vt (4.1) 

and outputs the potential 

u?:=5*A?. (4.2) 
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Using a simple Laplace transform argument and the estimates of [1] (see also [10]), it is 
easy to prove that (14. ip has at most one continuous causal solution. 

Before proceeding to state and prove the main result of this section, we are going to 
introduce some constants related to the geometry of the problem and associated functional 
inequalities. 

4.1 Some inequalities 

Let R > be such that 

07 C B := B(0; R) = {x G R d : |x| < R} 

and let us consider the balls Bt := B(0; R + T) for T > 0. Let then Ct > be taken so 
that 

\\u\\ Bt < C t \\S7u\\ Bt Vm G Hl(B T ). (4.3) 

A simple scaling argument shows that we can take Ct = Cq(1 + T/R). Therefore, the 
constant Ct grows linearly with T . The trace operator on the boundary OB? will be 
denoted jt- 

We will also consider a constant for the following two-sided trace inequality: 

IMIi/2,r < C r |M|i,B G H l (B ). (4.4) 

Next, we consider a one-sided lifting of the trace onto T in the form of a bounded linear 
operator L : H l / 2 (T) H\R d \ T) such that 

Ly? = in fi+ and ^ = <p Wip G F 1/2 (r). (4.5) 

Note that = V 9 an d that there exists Cl > such that 

||^||i iRd \ r = ||^||i,n_ < Cx|^||i A r G tf^r). (4.6) 

Finally, using the weak definition of the normal derivative, we can fix a constant C v > 
such that 

IIEMII-W < C,(||A M ||| fAr + ||V M ||| oVr ) 1/2 Vn G #A(£ \ T). (4.7) 

4.2 Estimates for the Galerkin projection 

Theorem 4.1. Let A G Wq(M; ff~ 1 / 2 (r)). T/ien £/ie solution of the semidiscrete problem 
04. ip and its associated potential (14. 2 p satisfy 

e C (M; ^ 1/2 (r)), u g g Cq(1R; tf 1 ^)). 

Moreover, for all t > 0, 

iK^Ux,^ < 2C r (||A(t)||_ 1/2ir + 7lT^5- 1/2 (A,t)), (4.8) 
||A£(*)||- 1/2j r < (l + C'rC' I/ )||A(t)||_i /2j r + V^C'raS2 1/a (A,0, (4-9) 

^ V2 (A,t) := jf (||A(r)||_ 1/2 , r + ||A(r)||_ 1/2 , r )dr. 
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The proof of Theorem 14.11 will occupy the remainder of this section. The proof will 
never use that Xh is finite dimensional. If we take Xh = i7 _1 ' 2 (r), ( 14. 8 p gives a bound 
for the single layer acoustic operator that reproves [8j Theorem 3.1]. 

Theorem |4]T] will be proved for A G Cq(M; if _1 / 2 (r)). The extension to the general case 
follows by a simple density argument. Also, we will prove the results in a finite interval 
[0,T], and bounds (j4.8p - (l4.9p will only be proved for t = T. Since T is arbitrary this is 
equivalent to having proved the results for any t. 

Because of the finite speed of propagation of solutions to the wave equation, it is 
possible to understand (formally at the beginning) as a solution of the following wave 
propagation problem on a truncated domain with non-standard transmission conditions 
(see USES]) 

v%(t) G Hq(B t ) < t < T, (4.10a) 
u%(t) = Au%(t) < t < T, (4.10b) 
juf(t) - (V * A)(f) eX° h < t < T, (4.10c) 
[cU«(i)] eX h < t < T, (4.10d) 
u£(0) = r£(0) = 0. (4.10e) 

The set X% is the polar set or annihilator of X^, i.e., 

X° h := {p G H l ' 2 {T) : (fi h ,p) r = V/x fc G X,}. 

If u G H\(Bt \ r), the condition \d v u\ G X^ can be rewritten as a set of restrictions 

([M.P>r = V P GX°, 

or equivalently, as 

(Vu, Vw) BT \r + (Am, v) BT \ r = Vf G V^ T (4.11) 

where 

:= {v G Hq(B t ) : (4.12) 

Proposition 4.2 (Uniqueness). Problem (I4.10p has at most one solution 

ut G C 2 ([0, r]; ^ 2 (5t)) n C x ([0, T]; H£(B T )) n C([0, T];Hi{B T \ Y)). (4.13) 

7/ i/izs solution exists and A^ := [(Vu^], £/ien and A^ coincide with the solution of 
dUD-dl^D on #ie mieruaZ [0,T]. 

Proof. Uniqueness of solution follows from an elementary energy argument. Careful, but 
not complicated, use of extension operators from B? to M. d and from [0, oo) to R in the 
sense of vector valued distributions, can be used following [SI Sections 4 & 5], to prove the 
relation between a strong solution of the transmission problem and a weak distributional 
solution of (P1)- (Q]) . □ 
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4.3 An underlying initial value problem 

Associated to the space V£ given in (I4.12p . we consider the space 

D T h := {v G V? : Av G L 2 {B T \ V), \d v v\ G X h } (4.14) 
= {ve V? n Hi(B T \ T) : (Vv, Vw) Bt + (Au, w) BT \r = G } 

(both definitions coincide by (14. lip ). In the frame of the triple C C L 2 (Bt) we 
can consider the unbounded operator A : D\ C L 2 (B T ) — >■ L 2 (B T ) (the distributional 
Laplacian in 5^ \ T)) and the initial value problem 

ii(t) = Av(t) + /(t) Vt > 0, t»(0) = u(0) = 0. (4.15) 

Strong solutions of this problem will be those with values in D\. Following P, Section 8 
and Appendix], we can also consider weak solutions of (I4.15p . In order to do that, we start 
in the Gelfand triple C L 2 (B T ) = L 2 (B T )' C (V^)', we extend A to an unbounded 
operator A : V h T C (V h T )' — > (V^)' and we only look for solutions with values in V£ . In 
this case, the additional boundary condition that appears in the definition of D\ and the 
extended definition of the Laplace operator are part of the same expression: namely, the 
differential equation (I4.15P is to be understood as 

(v(t), w) (V Ty xV T + (Vv(t), Vw) Bt = {f{t),w) BT Vw G V£. (4.16) 

The following two results follow from [U Appendix]. 

Proposition 4.3 (Strong solutions). If f : [0, oo) — >■ V£ is continuous, then there exists 
v : [0, oo) — > D\ such that 

v G C 2 ([0, oo); L 2 (B T )) n C\[0, oo); H^Bt)) n C([0, oo); H\{B T \ T)), (4.17) 

satisfying ( 14.15ft and the bounds 

\W)\\b t < C T f \\f(r)\\ BT dr, 
Jo 

\\Vv(t)\\Br < [ \\f(r)\\B T dr, 
Jo 

\\Mt)h T \r < f \\Vf(r)\\ BT dr. 
Jo 

Proposition 4.4 (Weak solutions). If f : [0, oo) — » L 2 (Bt) is continuous, then there 
exists v : [0, oo) — > such that 

v G C 2 ([0, oo); (V?)') n ^([0, oo); L 2 (B T )) n C([0, oo); H*(Br)), 

satisfying the weak form of ( 14.151) (see ( 14. 16ft ) and the bounds 

C^\\v{t)\\ BT <\\Vv{t)\\ BT < f ||/(r)|| Sr dr. 

Jo 

Moreover, v(t) := f^v(r)dr is a continuous D\ -valued function. 
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4.4 A decomposition of 

As explained in Section I4.2[ we are going to deal with A G C 2 ([0,T]; H' 1 / 2 ^)) such that 
A(0) = A(0) = and obtain bounds in the interval [0,T] as well. The proof of this result 
follows from the decomposition 

u% = S*\ + u° h + vl (4.18) 
where u° h : [0, oo) — >■ Vj[ solves the steady-state transmission problems (for all t > 0) 

u h (t)eV?, -Au° h (t)+u° h (t) = mB T \T, [d v u° h {t)]+\(t)eX h , (4.19) 
and v ° : [0, oo) — >• £)j[ is a solution of the evolution problem 

v° h (t) = A«g(f) + <(i) - t$(f) = At#(f) + Attg(t) - t#(f) t > 0, 

4.20 

«g(0) = tig(O) = 0. 

We start by analyzing the three terms in (14. 18ft one by one. Note that we still need to 
show that the decomposition ( 14.181) holds true, that is, that the sum of the three functions 
in the right hand side of (14.181) is u®. 

1. By Theorem 3.1 and Proposition 4.2 in [8], it follows that 

S * A G C 2 ([0, T}; L\B T )) n ^([0, T}; H]{B T )) n C([0, T];H^(B T \ T)), (4.21) 
that for all t G [0, T], 

£(S * A)(f) = A(S * A)(t), {d v (S * A)(f)] = \(t), 
(S*A)(0) = ±(S*A)(0) = 0, 

and 

||(5 * A)(0||i, Bt < C r (||A(0||-i/ 2j r + vT+C«B 2 - 1/2 (A, *)) < t < T. (4.23) 

2. We next analyze the behavior of u° h . The variational formulation of (I4.19P is 

u° h (t)eV?, (Vu° h (t),Vv) BT + (u° h (t),v) BT = -(X(t), 1 v) r VveV?. (4.24) 

This is a well posed problem, that depends on t, only because data depend on t. In 
particular, 

u° h G C 2 ([0, T];Hq(B t ) n H\{Bt \ T)), <(0) = <(0) = 0. (4.25) 
The variational formulation ( 14. 24ft and the trace inequality (14. 6 p show that 

\\<{t)h*r <^r||A(0 H-iAi, (4.26) 
On the other hand, since Au°(t) = then ( 14. 7p implies that 

HM(0]||-i/v < C v \\u° h (t)\\i,B T < aC r ||A(t)||_ lAr . (4.27) 
Differentiating ( 14.231) twice with respect to t, it also follows that 

||^(t)||i,B T <Cr||A(t)||_ 1/2 , r . (4.28) 
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3. We can apply Proposition 14.31 to the solution of (I4.20p . taking / := u\ — vP h . It then 
follows that 

v° h e C 2 ([0, oo); L\B T )) n C\[0, oo); H l {B T )) n C([0, oo); #i(B T \ T)). (4.29) 
We also obtain the bounds 

KWIIflr < a T K(r)-fi!l(r)|MT < C T C r B 2 1/2 (X,t) (4.30) 

(we have used (I4.26P and (I4.28P in the last step), 

\\Vv° h (t)\\ BT <C r B 2 1/2 (X,t) (4.31) 

and 

\\Av° h (t)\\ BT \ r <C r B 2 1/2 (\,t). (4.32) 

From the bound for the jump of the normal derivative (14.71) and (I4.31l) -( j4.32l) . it follows 
that 

IIMMHI-W < V2C r C u B 2 1/2 (X,t). (4.33) 
4.5 Proof of Theorem 14.11 

Let now be defined by (I4,18p , To prove that satisfies (14.131) . we just have to use 
(I4.2ip . (I4.25P and (I4.29p . To prove that satisfies problem (14.101) . we add the equations 
that are satisfied by the three components of the sum, namely (14.221) . (I4.19P and (I4.20p . 
The uniqueness result of Proposition 14.21 shows then that the function defined by (14.181) 
can be identified with the solution of (I4.ip - (I4.2I) in the interval [0,T]. 
By (14^31) . fl426|) . fl430|) and (jPTjl . it follows that 

||«? WIIi,b t < 2C r ||A(t)||_ 1/2 ,r + C r (>/l + C 2 t + yjl + C^B 2 1/2 (\,t). (4.34) 

Taking t = T and using that T is arbitrary, (I4.8p follows. 

Since \%(t) = = X(t) + [d v ul{t)\ + [d v v° h (t)] t and \d v -\ : H^(Bt \ T) 

//"-^(r) is bounded, then f l4T25j) and ( jOSD imply that A^ G C([0, T), i?-V 2 (r)). The 
uniqueness argument of Proposition 14.21 then proves that the solution of (14. ip is a causal 
continuous if ~ 1//2 (r)-valued function. Finally, inequalities (I4.27P and (I4.33P prove (14. 9p . 

5 The semidiscrete Galerkin solver 

In this section we study how Galerkin semidiscretization depends on data. Our starting 
point is a causal function (p :R—t H 1 / 2 ^). We then consider the function A^ : R X h 
such that 

(fi h , (V * A£)(t)) r = (fi h , \tp[t) - {K * <p)(t)) r Vfi h e X h (5.1) 
and the associated exterior solution 

u£ := -5 * A£ - Z> * ¥>. (5.2) 
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Using a simple Laplace transform argument and the estimates of [1] (see also [10]), it is 
easy to prove that (15.11) has at most one continuous causal solution. Moreover, uniqueness 
can be also established for weaker solutions, where for instance, is the distributional 
derivative of a continuous causal A^-valued function. 

5.1 Estimates for the Galerkin solver 

Theorem 5.1. Let cp E Wq(K; if 1//2 (r)). Then the solution of the semidiscrete problem 
(15. ip and its associated potential (15 .2p satisfy 

\t e C(R; #~ 1/2 (r)), ul e C X (M; H 1 (M d \ r)). 

Moreover, for all t > 0, 

K(*)IW < C L (\\ V (t)\\y 2 , r + VT+Cl 5 2 1/2 (^,t)), (5.3) 
K(t)||_ 1/2 , r < V2C v C L (4\\<p(t)\\i/2,r + 2\\(p{t)\\ 1/2jT + B\ /2 {y,t)), (5.4) 

Bl /2 {^t) := ^(||^(r)|| 1/2)r + ||^(r)|| 1/2 , r )dr, 

S 4 1/2 (^,t) := 4 J B 2 1/2 (^,t) + J B 2 1/2 (^,t). 

The proof of this result follows partially the steps of the proof of Theorem 14.11 The 
analysis will be more involved because of the occurrence of a non-homogeneous essential 
transmission condition (see ( I5.5dl) below), that cannot be easily lifted with a continuous 
potential. Like in Section HJ we will prove the estimates on a fixed time interval [0, T], 
taking advantage of the fact that finite speed of propagation will allow us to impose a 
homogeneous Dirichlet boundary condition in dB^. Again, we will only assume that 
(p E C A ([0,T\;H^ 2 (T)) with p( fc )(0) = for k < 3. The result for general <p can be 
extended with a density argument. 

One of the keys towards the proof of the result lies in the fact if is smooth enough, 
then is a strong solution of the following problem: 



ut(t)EH\B T \T) 


< t < T, 


(5.5a) 


lTut{t) = 


0<t<T, 


(5.5b) 


Kit) = Au?(t) 


< t < T, 


(5.5c) 


[X(t)] = ^) 


< t < T, 


(5.5d) 




< t < T, 


(5.5e) 


[d v ut{t)\ e x h 


< t < T, 


(5.5f) 


<(0) = <(0) = 0. 




(5.5g) 



At the same time, we can combine (I5.5cl) and ( I5.5fl) . multiply by w G and integrate, 
to obtain a weaker form of the differential equation and the natural boundary condition 

{u%{t),w) BT + {Vul(t),Vw) BT \ T = V«;6l4 T , < t < T. (5.6) 
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If we now define 

ft 



«t) ■= / <(r)dr, (5.7) 

Jo 

it also follows that 

K(t), W ) BT + (V<(t),V W ) BT \ r = VweV?, 0<t<T. (5.8) 

The following double uniqueness result will help us recognize in the two decompositions 
that will be given below. 

Proposition 5.2 (Uniqueness). There exists at most one 

ut G C\[0, T}; L\B T )) n C([0, T]\H l (B T \ V)) 



satisfying the essential boundary and transmission conditions (15.5b|) . (15.5d|) . and (I5.5el) . 



the initial conditions ( |5.5g[ ) , and snc/i that w% defined with (15.71) satisfies (15. 8p . If such 
a solution exists and, in addition, w% G C([0, T]; H\(B T \ Y), then Aj£ = — 3t[c?i/w£] (with 
differentiation in the sense of H -1 ' 2 (T) -valued distributions) solves (15. ip and satisfies 

(E2D. 



Proof. Let be a solution of ( I5.5bl) .( l5"".5dl) with (p = 0, (I5.5e[) and ( |5.5gP such that the 



corresponding w\ satisfies ( 15. 8p . Then : [0,T] — >■ V h T is continuous and we have 
enough regularity to test ( 15. 8p with u^~°(t) = w^~°(t) G and prove that 

^(lK =0 WIII T + l|v< =0 (t)||| TVr )=o. 

The proof of uniqueness of solution is now straightforward. The final statement follows 
from the same arguments developed in 0, Section 5]. □ 

Corollary 5.3. There exists at most one 

ut G C\[0, T);L 2 (B T )) n C\[0, T]\ H\B T \ V)) n C([0, T]; /Ti(B T \ T)) 

t/iat solves ( 15. 5ft . Moreover, if such a solution exists, then := — [c^ttjO solves ( 15. ip and 
satisfies ( 15.2ft . 

5.2 A first decomposition of 

For the arguments of this section, we only need </? G C 2 ([0, T]; i? 1//2 (r)) with y>(0) = 
0(0) = 0. In a first step, we formally decompose 

< = < + < (5.9) 
where u° h : [0, T] — > H 1 (Bt \ T) is the solution of the variational problems (for < t < T) 
u° h (t) E Hi(B T \T), 

lT u° h (t) = 0, [ 7 <(t)] = <p(t), i + ul(t) G XI (5.10) 
(Vu°(f), Vt^yr + «(*), =0 Vf G 
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and v® : [0, T] — > is a solution of 

KWMtyfyxv? + (V^(t), Vw) Bt 

= (u° h {t)-u° h (t),w) BT VweV^ 0<t<T, ( 5 . n) 

v° h (0) = «g(0) = 0. 

1. Problem (I5.10P has a unique solution by a simple coercivity argument. Since depen- 
dence on t happens only through the non-homogeneous essential transmission condition 
(compare with (14. 24ft . where the condition was natural), it is simple to prove that 

u° h £ C 2 ([0, T]; H\B T \ T)), „°(0) = u° h (0) = 0. (5.12) 

Using the lifting operator L given in (|4.5p -( PO)]) and taking u Q h {t) — Lip{t) £ as a test 
function in ( I5.10p . we can prove that 

lK(t)lkBr\r<CL||v(t)|| 1 / ai r. (5.13) 

Taking second derivatives with respect to time in (15.101) and the using same kind of 
argument, we prove 

||4Wlli,^\r<^L||¥3(t)||i/ 2 ,r. (5.14) 

2. Proposition 14.41 can now be invoked to prove that the evolution problem ( 15.111) has a 
unique solution 

v° h £ C^O, T]; Hl(B T )) n C([0, T]; L 2 (B T )), (5.15) 

satisfying the bounds 

C T l \K{t)\\B T < \\Vvl{t)\\ BT < f\\ul{r)-ul{r)\\ BT dT<C L Bl /2 {^t), (5.16) 

Jo 

where we have used (I5.13I) - (I5.14I) in the last inequality. 
5.3 A second decomposition of 

An alternative decomposition to (15. 9p is needed to bound the density A^. From this 
moment on, we need cp £ Cq([0,T]; i7 1 / 2 (r)). We now write 

ut = u) l + vl (5.17) 

where u\ : [0, T] — > H l {Bx \ T) is the solution of the variational problems (for < t < T) 

ul(t)EH\B T \T), 

lTul(t) = 0, bu\{t)] = ip{t), 7 + <WeX°, (5.18) 
{Vu\{t), Vv) BAr + {ul{t),v) BT = {L(p{t) - Ltp{t),v) BT \/v £ Vf , 
and v\ : [0, T] — > D\ is a solution of 

vl(t) = Avl(t) + f(t), 0<t<T, 

5.19 

vi(o) = m = o, 

where / := u\ — u\ + Lip — Lip = Au\ — u\ : [0, T] — > is continuous. 
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1. It is clear that u\ G C 2 ([0, T\\ H\{B T \ V)) and ^(0) = ii l h (0) = 0. Also, using 
v = u\(t) — L<f(t) G as test function in (15.181) . we can prove that 

K(0lli,B T \r < ||M0lli,B T \r+ 11^(0 - L<p{t)\\ ltBT \r (5.20) 
< C L (3||^(t)|| 1/2 ,r + ||v?(t)||i/2,r) 

and therefore 

||A4(t)|| BT \ r <C L (4|^(t)|| lA r + 2||^(t)|| lA r) (5.21) 

2. Regularity of the solution of (15. 19[) is given by (I4.17P in Proposition 14.31 Since (15.201) 
can be differentiated twice with respect to time, we can bound 

||V/(t)|| BT < K(t)||i >BTV r + H(t)\\iA-\v + C L (\\<p(t)\\w + II0(*)II W) 
< C L (4||^(t)|| lAr + 5|| ¥ 3(t)|| 1/2 , r + ||^ 4 )(t)|| lAr ) 

and thus 

||A^|| BAr < / \\Vf{r)\\ BT dr < C L B 4 { V ,t) < t < T. (5.22) 
Jo 

5.4 Proof of Theorem 15.11 

We first define v% with (I5.9p . where u° h solves (15.101) and t>° solves (15.111) . By (15.101) . 
(15. lip . (I5.12p and (15.151) . it follows that uf t is the only weak solution of (15. 5p in the sense 
of Proposition 15.21 and it can be thus identified with the solution of (15. ip ~ ( 15721) . As a direct 
consequence of (15.131) and (15.161) . it follows that 

\K(t)\\ 1>BT \T < C L (\\<p{t)\\ 1/2) r + ^/l + C|,S 2 1/2 (^t)), < t < T, (5.23) 

and 

||V<(t)|| BT < C L (\\<p(t)\\ w + 5 2 1/2 (^,t)), < t < T. (5.24) 

We next define vif h = u\ + v\, where u\ satisfies (15 . 181) and v\ satisfies (I5.19p . It is 
clear that is a strong solution of (14. 3p and therefore (Corollary 15.31) coincides with 
the solution of dEID-Q. Since u£ G C([0, T}; H^(B T \ T), then A£ = -[d„v%\ G 
C([0,T];i/- 1 / 2 (r)). By (OOl and (OT]l it follows that 

\\Auf L (t)\\B T \r < C L (q<p(t)\\w + 2\\cp(t)\\ 1/2jT + 5 4 1/2 (^,t)) < t < T. (5.25) 

This inequality, (15.241) . and (14 .7p prove finally that 

||A£(t)|U A r < VC L C u (4\\<p{t)\\ 1/2>r + 2|Ki)|| 1/2>r + B^^t)) 0<t<T. 
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5.5 Proof of Theorem 13.11 

We use the notation (14. ip for the Galerkin projection and (15.11) for the Galerkin solver. 
Since, (LT^A)^ = IT^A, we can bound 



|| (A - X h )(t)U /2 ,r < ||(A - n,A)(t)||_ 1/2 , r + ||(A - n fc A)£(t)||_ 1/2 ,r + \\\™ 
The bound (13.51) follows then from Theorems 14.11 and 15.11 Similarly we write 



-i/2,r- 



<f>— <Ph 
h 



\\{u - u h )(t)\\ hRd < \\S * (A - n h A)(*)|| ljR d + ||5 * (A - n A A)^(t)|| ljR d + ||u 
and use 0, Theorem 3.1], Theorem 14.11 and Theorem 15. II to prove (I3.6p . 

6 Analysis of time discretization 
6.1 A passage to the Laplace domain 

For any value s G C + := {s G C : Res > 0}, we consider the fundamental solution of the 
differential operator A — s 2 , namely, 

f itf^Hx-yl), ford = 2, 
$(x,y; S ) := 1 e~ s ^ ior d = 3 

I 47r|x — y|' 

(The function Hq is the Hankel function of the first kind and order zero.) We also 
consider the single and double layer potentials 



S(s)X:=J^( ■ ,y; S )A(y)dr(y) : K \ T — >• C, 

D(s)<p := J d v(y) $( ■ ,y;sMy)dr(y) : R d \ T -> C 

and the associated integral operators 

V(s) := 7 + S(s) = ^S(s) K(s) := ^ + B(s) + § 7 ~D(s)- 

Consider then the Laplace transforms of the semidiscrete data $^ := £{<fh} and of the 
semidiscrete solutions := C{\h} and := C{uh}- For z G C + and G G H 1 ^ 2 ^), we 
consider the uniquely solvable transmission problem looking for V G ff 1 (M d \r) such that 

z 2 V-AV = inM d \r (6.1a) 

[ 7 V] = G, (6.1b) 

fV6l2, (6.1c) 

[d„V\ G X h . (6.1d) 

Proposition 6.1. For all s G C +; U^(s) is the unique solution of the transmission prob- 
lem (16.11) with z = s and G = $/i(s). Moreover A^(s) = — \d v U^(s)] . 
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Proof. Note that for all s G C+, Ah(s) G Xh and that we have the relationships 

V(s)A h (s) - + K(s)$ ft (s) G U A (s) = -S(s)A fc (s) - D(s)$ h (s). (6.2) 

This proves the result. □ 



The CQ discretization affects all four convolutions in ( 16. 2ft . It defines causal functions 
ul and A£ such that u%(t n ) = u%[ri\ and A^(t n ) = A£[n]. Let then A£ := £{A£} and 
U£ :=£{<}. 

Proposition 6.2. For all s G C +; U£(s) is the unique solution of the transmission prob- 
lem (16.11) with 

5(e- SK ) 

z = s K := G C + Vs G C+, 

and G = $ h (s). Moreover A£(s) = -[0„U£(s)]. 
Proof. By construction A£(s) G and 

V(s K )A£( S ) - + K(s«)$ h (s) G X°, U£(s) = -S(s K )A K h {s) - D( Sli )<S> h (s). (6.3) 

The result is then straightforward. □ 

Note that each occurence of s K in ( 16. 3 j) corresponds to the discretization of a convolu- 
tion process. 

We finally consider the errors of time discretization of the semidiscrete-in-space prob- 
lem 

e = el := u h -< E := U„ - U£ = C{e} e = e K h := A, - AJJ. 
On time steps, we will be considering the errors 

e„ := e(t n ) = u h (t n ) - u%[n], e n := e(t n ) = X h (t n ) - A£[n]. 
Applying Propositions 16.11 and 16.21 it follows that for all s G C + 

E(s) G D h := {u G H\R d ) fl tfi(M d \ T) : 7 u e X°, [0„u] G X h }, (6.4) 

and 

#3( a ) - AE( S ) = (si - , 2 )U,( S ) =: Q(s). (6.5) 

In the time domain, the function 8 = 8%, whose Laplace transform is 0, corresponds to a 
consistency error of the time discretization -the approximation of the second derivative by 
the particular CQ scheme-, applied to the semidiscrete-in-space solution. Before we start 
the analysis of each of the time discretization methods, let us mention the integration-by- 
parts formula in Dh, which will play an important role in the forthcoming analysis: 

(Au,v) K d\ r + (Vu, Vv) K d = {{d v u\,~jv) = Mu,veD h . (6.6) 

Bounds with respect to data will be given in terms of the following quantities: 

^ 7 W):=E / hhWhftfdT (*>5). 
The following product (semi) norm 

\\(u,v)\l := (llVnll^ + I^H^' 172 
will be used to simplify some formulas. 
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6.2 Analysis of the Backward Euler discretization 
Proposition 6.3. For all n > 1, 

IKe^i^-e^))! < KCt n (l + t n )Bl /2 (ip h ,t n ). (6.7) 

Proof. Let f n := -(e n — e„_i). Noticing that for the Backward Euler discretization 
s K = -(1 — e~ SK ), the error equation (16.51) can be written as 

e n - e n _i = «/„, /„ - / n _i = /tAe n + k0„ (6.8) 

where we can bound the consistency error as 



u h (t n ) - 2u h (t n - 1 ) + Uh(t n - 2 ) _••/.> 

2 "H"n > 



< ^ max ||4 3) (*)|| Rd . (6.9) 

O t n -2<t<t n 



Testing the equations (16. 8ft with — Ae„ and / n respectively, adding the results, and ap- 
plying the integration by parts formula (16. 6ft (note that e n G -D/i for all n, since e takes 
values in this space by (16.4jl ). it follows that 

(V )]R<* + (/n, /n)lRd = (Ve 

)]R d + (/nj /n-l)lR d + K (@n, fn)n. d 

and therefore 

II (en, /«) II < lll(en-i, /n-i)||| + « Vn. (6.10) 

Then, by induction 

n 

|||(e n ,/n)||| < ll^'llK d - 

Using now (16. 9p and Theorem 15.11 (recall that (p i— >■ is a convolution operator and 
therefore commutes with differentiation), it follows that 



< ^^(^ ||ri 3) (t)||i/ 2 ,r + ^/l + CgB 2 1/2 (ri 3) ,t n )), (6.11) 



where we have also used the fact that the constant Q of (14. 31) grows with t. Adding the 
bounds (16. lip for different values of j, using the overestimate 

max ||ri 3) (t)|| lAr < max ||ri 3) (t)|| 1/2 ,r < / K\t) ||i A rdt, 

the result follows. □ 
Theorem 6.4. For all n > 1 



|e w ||i,Md < Kt n (l + t n )B s ' (tph,t n ) (6.12) 

-i/2,: 



Ml-i/ 2 ,r < K(l+t 2 n )Bl /2 (<p h ,t n ). (6.13) 
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Proof. Using the notation of the proof of Proposition I6.3[ and since e n = e„_i + K,f n , we 
can easily bound (using Proposition 16. 3j) 

n 

\\e n \\Rd < H//'IU d - K( 3tn (1 + t n ) Bl /2 (if h ,t n ). 
i=i 

This inequality and Proposition 16.31 prove (I6.12p . 

Note now that f n = f(t n ), where / = (e — e( • — k))/k. Therefore, using the second 
of the equalities (16. 8p . we can bound 

< max ||/(t)|| R d + ||0 n || R a. (6.14) 

For t G t n ], we can construct a mesh with time-step k that includes t. Then, applying 
Proposition 16.11 with data (fih on this mesh, it follows that 

||/(t)|| R d < nCt(l+t) Bl /2 (ip h ,t) < KCt n (l + t n )Bl /2 (cp h ,t n ). 

This inequality, the bound (16.111) for the consistency error and (16.141) provide a bound for 
the Laplacian of the error 

|| Ae n || Rd \ r < kC (1 + t 2 n ) 5 6 1/2 (^, t n ). (6.15) 

Since by Propositions 16.11 and 16.21 it follows that [cte n ] = — e n , the bound (16 . 131) is a 
direct consequence of (14.71) . Proposition 16.31 and (16.161) . □ 

6.3 Analysis of the BDF2 discretization 
Lemma 6.5. Given G G L 2 (R d ), let V G D h solve 

s 2 K V- AV= (s 2 K -s 2 )s- 3 G inR d \T 

Then 

|||(V,s K V)||| <Ck 2 . \ sl ||G|| Rd VsgC + , Vk, V/i. 
mmjl, tie sj 

Proof. Consider first a general z G C+, H G L 2 (M d ) and V G -D?, such that z 2 V - AV = H. 
By ([61D, it follows that 

^ 2 ||V||^+||VV||^ = (H,V) Rd 

and therefore 

(Rez)||(V,zV)| < ||H|| Rd . (6.16) 
Taking now H := (s 2 , — s 2 )s~ 3 G, z = s K , and noticing that 

\s K \ < Ci\s\, \s K — s\ < C 2 k 2 |s| 3 , Res K > C 3 minjl, Re s}, Vs G C+,Vk, 

the result follows from (16. 16ft . □ 



|Ae n || R d\ r 



K 



■{fn — fn-l) — 9r, 
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In the next results we will refer to the operator 

d K g := i(§<7 - 2g(- - K ) + §<?(• - 2k)) C{d K g} = s K G(s), 
which is the discrete derivative associated to the BDF2 method. 

Proposition 6.6. Let f := d K e and f n := f(t n ) = -(§e n — 2e n _i + |e n _ 2 ). Then for all 
n > 1 

\\{e n J n )\\<K 2 Ct n {l + t 2 n )Bl /2 { Vh ,t n ). 
Proof. By (16. 5ft and Lemma 16.51 it follows that 

|||(E(5),F( S ))||| < Ck 2 \ S] \\s 5 XJ h (s)\\ Rd s e C + . 
minjl, tie s} 

Using then [HI Theorem 7.1], it follows that 

H(e(f),/(f))l < f H^WMt Vt. (6.17) 

Since by Theorem 15. 1[ we can bound 



|w h (r)|| R ddr < C(l +t 2 )5 2 1/2 (^,t), (6.18) 



the result is a direct consequence of (j6.17|) . □ 
Theorem 6.7. For all n>\, 

\\e n \\i^ < K 2 t n (l + t 3 n )Bl /2 (ip h ,t n ), (6.19) 

1/2 



kn||-i/ 2) r < /€ 2 C(l + 4)S 9 1/2 (^,t n )- (6-20) 



Proof. The proof is very similar to that of Theorem 16 .41 Using a simple stability argument 
for recurrences, we first show that 

n 

||e n ||»d < \\fj\W d - 
i=i 

This inequality and Proposition 16.61 prove (I6.19p . We next use that Ae n = (d K f)(t n ) — 6 n 
(see (16.51) and the definition of / in Proposition 16. 6p to bound 

||Ae„|| R d\ r < §||£(/ n ~ fn-l)\Wd + §|| £(/»-! - fn-2)\\Rd + \\(6%U h - U h )(t n )\\ Rd 



^ 2 + \\f( T )h* + CK max \\u\\t) 



Using then (16.171) and (I6.18P we obtain a bound 

||Ae n || Md \ r < k 2 (1 +t 2 n )Bl /2 ((p h ,t n ) 
and (I6.20p follows from this and Proposition 16.61 using that e n = — ld u e n }. □ 
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6.4 Analysis of the trapezoidal rule discretization 
Proposition 6.8. For all n > 1 

\l(l(e n+1 + e n ),±(e n+1 -e n ))\l < k 2 C t n+1 Bl /2 ((pf\ t n+1 ). 
Proof. In the case of the trapezoidal rule, the error equation (16. 5p can be written as 

4,(1 - e" SK ) 2 E( S ) = |(1 + e" SK ) 2 AE( S ) + \k 2 {\ - e— ) 2 U h (s) - ±(1 + e"« )VU h («). 
In the time domain, this gives 

4r (e n+ i - 2e n + e n _i) = jA(e n+ i + 2e n + e n _i) + x«, (6.21) 

where 

Xn ■= 4r(«fc(Wi) - 2u ft (t n ) + «fc(*n-l)) - + 2« ft (t„) + itfc(*r*— x))- 

Using the integration by parts formula (16. 6p we obtain 

^(e n+ i - e n ) - i(e n - e n _i), uj ^ (6.22) 

+ f (V(|(e n+1 + e n ) + V(|(e n + e n _i), Vv) = «(*„, i>)m" Vv G D h . 

Testing (16T22|) with 

v := ^(e n+ i - e n _i) = \{e n+ i - e n ) + ~(e n - e„_i) = ~ (|(e n+ i + e n ) - |(e n + e n _i) J , 
it follows that 

|||(|(e„+i + e n ), i(e n+ i - e„))|| 2 = |||(|(e n + e n _ x ), \{e n - e n _i))||| 2 

+ K (Xn, -(e n +i — e n ) + ~(e n — e n _i)) Rd 

n 

= *( e i+! ~ e i) + K e i ~ e i-i))iR d - ( 6 - 23 ) 

j=0 

Consider now the mesh-grid with nodes in the midpoints £- + i := (j + |)/t, the piecewise 
constant function \ such that = Xj in (tj-±,tj + ±), and the continuous piecewise 
linear functions / and e such that 

f(t j+ i) = \{f j+ i - fj), e(t j+ x) = i(e i+ i + e,-). 

We can then write (I6.23P as 

lll(^ + l),7(i n+ l))||| 2 = 2 [ tn+ \x(T)J(T)) Rd dT. 

2 2 Jo 

Given n, we choose n* < n such that |||(e(t n * + i ), f{t n * + i))f is maximized. Then 

lll(e(W§) J(Wi))I 2 < 2[(e(Wi)./(Wi))l Hx(r)||»-dT 
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and therefore (after comparing with the maximum at i n * + 1 and overestimating the integral 
in the right-hand side) 

l\(e(t n+ i)J(t n+ i))\i <2 [ n+i ||x(r)|| R ddr <2t n+1 max\\ Xj \U* Vn. 
2 2 Jo o<n 

Since 

llXnlkd < Ck 2 max ||4 4) (r)|| R d, (6.24) 
the result follows by (16.181) . □ 
Theorem 6.9. For all n > 1, 

\\e n \\ Rd < K 2 Ct 2 n B l 7 /2 {<p h ,t n ), (6.25) 
||i(e rH .i + 2e ft + £„_!) ||_ 1/2ir < k 2 C(1 + t 2 n+1 )Bl /2 (<p h ,t n+1 ). (6.26) 

Proof. Since 

n— 1 j 

e„ = KV-( 

the first bound follows from Proposition 16.81 By (I6.2ip we can bound 

|||A(e n+ i + 2e n + e n _i)|| R< i\ r < C max ||e(r)|| E d + ||Xnlk«i ( 6 -27) 

in-l<T<tn + l 

The first term in the right-hand-side of (I6.27P can be bounded using (16.251) applied to if>h 
using a time-grid with time-step k that includes the point where the maximum is attained. 
The second term of (I6.27P is bounded using (I6.24j) . The result follows from the fact that 
£ n = -ld„e n ]. □ 



n-l 

1 

{e j+l -e jh 



7 Final comments 

In this paper we have given a full analysis of the discretization with Galerkin is space 
and three particular instances of Convolution Quadrature in time of a direct formulation 
for the exterior Dirichlet problem for the wave equation. The full error estimates are 
the result of Theorem 13.11 for the semidiscretization in space process and Theorems 16.41 
(backward Euler). ISTTl (BDF2). and 16.91 (trapezoidal rule) for time discretization. 
An indirect formulation, i.e., 

V*£ = <y9 u = S * £ 

follows from very similar arguments. The Galerkin projection is the same and therefore, 
the analysis of the semidiscrete in psace system is a particular case of the results in this 
paper. The Galerkin solver (see Section [5]) is however slightly different. While its analysis 
is not needed for the semidiscretization in space, it is needed for the time discretization. 
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This analysis is likely to be extremely similar to the one given here. In terms of its Laplace 
transform, the semidiscrete problem is 

s 2 XJ h (s) - AXJ h (s) = in R d \ T 
[ 7 U h (a)] = 0, 
7 U h (a) - Hs) e XI 
[d v \J h (s)j G X h . 

with S^(s) = [9 I/ U/ l (s)J. This is a very similar problem (same kind of transmission 
conditions) to problem (16 .11) (see also Proposition I6.ip . In particular, the error equations 
to compare semidiscrete and fully discrete solutions (I6.4l) - (l6.5p are the same as in the case 
of the direct formulation and all the arguments of Section [H] hold, contingent to having 
proved the estimates of Theorem 15.11 adapted to the new kind of Galerkin solver. 

All the arguments that have been used in this paper can be easily extended to the 
case of the single layer potential for the elastic wave equation in any dimension. 

Much of the analysis of SectionsHJl6]can be done using estimates in the Laplace domain. 
That gives a more streamlined way of proving estimates, although they come with either 
worse constants (for growth in time) or with higher continuity requirements: see [El Section 
7] for a comparison of Laplace domain and time domain techniques applied to estimating 
layer potentials and integral operators. The analysis of semidiscretization in space using 
the Laplace domain can be adapted from the techniques developed in [10] . Analysis of 
convolution quadrature can then be carried out using the very general results of Lubich 
[T2"] applied to the semidiscrete operators. It has to be noted, though, that the analysis 
in [12] does not cover the case of the trapezoidal rule (the reference [3] circunvents this 
difficulty nevertheless), while the relatively traditional time-domain analysis of Section 
16.41 -based on understading the semidiscrete equations as a transmission problem and, in 
particular, on the integration by parts formula (16. 6ft -. is applicable. Similarly, the use of 
multistage convolution quadrature [I], [5] and variable-step convolution quadrature [IT] 
can be applied using estimates in the Laplace domain and it remains to be seen whether 
a time-domain analysis is practicable and produces different or improved results. 
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